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Abstract 

We consider a pair of one-parameter (a) families of generalized two-qubit determinantal Hilbert- 
Schmidt probability distributions, p a (\p PT \) and g a (|p|), where p is a 4 x 4 density matrix, p PT , its 
partial transpose, with \p PT \ G jgg] and \p\ S [0, ^g]. The Dyson-index-like (random matrix) 

parameter a is ^ f° r the 9-dimensional generic two-rebit systems, 1 for the 15-dimensional generic 
two-qubit systems,. . .. Numerical (moment-based probability-distribution-reconstruction) analyses 
suggest the conjecture that the Fisher information-a measure over a— is identical for the two distinct 
families. Further, we study the values of p a (0), the probability densities at the separability- 
entanglement boundary, with evidence strongly indicating that P2(0) = ^|f^ = 3 2 5 i7 n > P3(0) = 

lIP = 24 3 1 23 17 ' -P 4 (°) = 1 266i04 5 = ^~W^§f '' ' ' " Des P ite extensive results of such a nature, we 
have not yet succeeded-in contrast with the corresponding rational- valued separability probabilities 
flhttp://arxiv.org/abs/1301.6617D -in generating an underlying, explanatory ("concise") formula for 



p a (0), even though the corresponding denominators of both sets of rational- valued results are closely 
related, having almost identical (small) prime factors. The first derivatives p' a (0) are positive for 
ol=\ and 1, but negative for a > 1. 
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I. INTRODUCTION 



The nature of the boundary between separable and entangled quantum states is a subject 
of interest in the study of the "geometry of entanglement sudden death" [U Fig. l(ii)]. Our 
analyses will provide information concerning the density of states in the vicinity of this 
boundary. Additionally, they suggest a conjecture pertaining to the Fisher information of 
two pertinent one-parameter families of determinantal probability distributions. Let us now 
place our study, and the accompanying results, in the context from which they emerged. 

Zyczkowski, Horodecki, Sanpera and Lewenstein, in a much cited 1998 article [2], posed 
"the question of how many entangled or, respectively, separable states there are in the set 
of all quantum states". In [3], we, in essence, provided a "concise" and highly generalized 
answer-though one still lacking a formal, rigorous demonstation-to this fundamental ques- 
tion in the context of two-qubit systems endowed with Hilbert-Schmidt (Euclidean/flat) 
measure [U [5]. The formula-for the Hilbert-Schmidt probability P(a) that a generic 2x2 
quantum system is separable-that was reported there took the form 

P(a) = + (1) 



where 



f(a) = P(a) 



p(a | , r(a)2-^-«r(3a+|)r(5a + 2) 
1 + ' 3r(a + l)r(2a + 3)r(5a+f)' K) 



and 



r(a) = 185000a 5 + 779750a 4 + 1289125a 3 + 1042015a 2 + 410694a + 63000 = (3) 
at 5a(25a(2a(740a + 3119) + 10313) +208403) +410694] +63000. 



Here a is a Dyson-index-like (random matrix) parameter that takes the value | for the 
9-dimensional two-rebit systems, 1 for the 15-dimensional two-qubit systems, and 2 for the 
27-dimensional two-quater(nionic)bit systems,. . . . (We have the relation n = 3 + 12a, where 
n is the dimensionality of the associated generic generalized 4x4 density matrices.) 

These separability probabilities were computed-in accordance with the seminal Peres- 
Horodecki results [HlEj-as the cumulative probabilities over the nonnegative interval |p PT | G 
[0, 2§g] of probability distributions p a (\p PT \), with \p PT \ G [— ^, ^g]. For a = |, 1 and 2, 
the formula yields (to arbitrarily high precision) P(a) — || « 0.453125, J| 0.242424 (cf. 
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FIG. 1: The a-parameterized family of probability distributions p a {\p PT \) over [—js, ^g] estimated 



on the basis of the first 245 Hilbert-Schmidt a-specific moments of \p |, given by ([4 

P eq. (B7)] sec. VII]) and ^ « 0.0804954, respectively. In the calculations of these 

quantities, the moment formula [PUl sec. D.4] 

prinv n!(a + l) n (2a + l) n 
^ l> 2an(3a + |) re (6« + |) 2n + U 

(_2n-l-5a) n (a) B (a+|) n f / -2=1, -n, a + 1, 2a + 1 \ 

2 4 ™ (3a+ |) n (6a+ |) 2n 5 4 \1 - n,n + 2 + 5a, 1 - n - a, § - n - a J 
was employed in a (Mathematica-implemented) Legendre-polynomial-based probability- 
distribution reconstruction procedure of Provost [11] (cf. |12j). 

Other than the separability probabilities given by the concise formula 0-((3]) for each 
value of a, however, nothing additional appears to be specifically known about the a- 
parameterized family of probability distributions over [— ^g], which we denote as 
Pa(\p PT \)- (Fig. [l] is a plot of this family, estimated on the basis of the first 245 Hilbert- 
Schmidt moments Q of |p PT |.) 

We have, thus, engaged in a research program intended to increase our knowledge of these 
probability distributions, which are of clear major quantum- information-theoretic interest. 
Our long-range goal is to fully characterize them-that is, develop explicit formulas, if pos- 
sible. From such formulas, one should then be able to derive/confirm the already-reported 
separability probabilities [3] as the cumulative probabilities of p a (\p PT \) over |p PT | G [0, ^]. 

Let us observe that twice the fourth root of the absolute value of the determinant |p PT | 
is a two-qubit "monotone under pure local operations preserving dimensions and classical 
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FIG. 2: Logarithms of y-intercepts p a (0) of estimated (based on 3,000 moments ([4])) Hilbert- 

16' 256 J 



Schmidt generalized two-qubit probability distributions over \p PT \ £ r 1 J 



communication and provides tight upper and lower bounds on concurrence" [I3j ec L- (5)]- 
(If |p PT | > 0, this "determinant-based measure" is taken to be zero.) 



PT\ 



II. PROBABILITY DENSITY ESTIMATES AT \p 

Possibly the most natural additional question to be first addressed is what is the value of 
the "y- intercept" of the probability distribution-that is the value the curve/function attains 
at |p PT | = (the separability-entanglement boundary). Such analyses would appear to be 
of relevance-among other areas-in studying the geometry of entanglement sudden death 

To begin, we estimated the y-intercepts jointly for the seventy integral and half-integral 
values a = \, 1, §, . . . 35, using the corresponding first 3,000 moments of |p PT | PH sec. D.4]- 
given by (|4|) in the Legendre-polynomial probability-distribution reconstruction procedure 
[TT] . The logs of these seventy estimates are displayed in Fig. [2j (A least-squares fit to this 
highly linear plot is provided by 7.35633 — 0.823151a.) 

Further, in Fig. [3j we show the analogous plot (exhibiting interesting non-monotonic 
behavior), based now on the estimated probability distribution of |p||p PT | over the interval 

[ _ H(i592> 2^62]~ ratner tllan \P PT \ OVer [ _ l6' ski" ( We n ° te tllat 110^92 = 2 ~ U ' 3 ~ 3 ') Hele 

2,600 moments pH sec. D.4], 

/, m I pt|"\ ( 2n ) ! ( 1 + «) 2 n( 1 + 2 «) 2 » f ~n,a,a+ |,-4n- 1 -5a \ 
{lPl lP ' ^ 2^ (3a + l) 2n (6a + |) 4n 4 3 V -2n - a, -2n - 2a, \ - n ' J ' 
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FIG. 3: Logarithms of y-intercepts of estimated (based on 2,600 moments ([5])) Hilbert-Schmidt 
generalized two-qubit probability distributions over |p||/0 PT | £ [— nsW' 2bW1 

were utilized in the reconstruction procedure. Rates of convergence are much slower in this 
"balanced" moment case than in the previous "unbalanced" case, so our further analyses 
here will fully focus on the unbalanced instance. 

In Table [j] we give some of our early estimates for the ^-intercepts (p a (0)) as a function of 
a-and, in some cases, advance apparently exact rational values for them. We are strongly 
confident-based on still further systematic computations with greater than 10,000 moments- 
in knowing the exact values for a = 2, . . . , 46. Nevertheless, Mathematica has not yielded 
an explanatory formula-as it accomplished in the separability-probability analysis with a = 
1, . . . , 28. Interestingly, most of the prime factors of the denominators of these exact results 
are the same as the prime factors of the already-reported corresponding exact separability 
probabilities. By way of illustration, the denominator of pi2(0) (Table [T]) is 3 2 ■ 79 • 83 ■ 89 ■ 
97 • 101 • 103 • 107 • 109 • 113, while that of the separability probability for a = 12 is identical 
to this, but for the replacement of 3 2 by 3 ■ ll 2 . The largest prime factors of p a (0) never 
exceed those of the corresponding separability probabilities, and can, on occasion, be less. 

Rates of convergence markedly increase as a itself increases, being weakest, unfortu- 
nately, in two (two-rebit and two-qubit) instances of immediate interest, a = | and 1. 
Also, estimates converge monotonically upward for a < 10, as the number of moments in- 
putted is increased, so, it is clear, that the corresponding reported results provide lower 
bounds on true values. (It has been suggested that the use of a Jabobi-polynomial-based 
probability-distribution reconstruction procedure [HI p. 742] might exhibit superior conver- 
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FIG. 4: Estimates (relatively slowly converging)-as a function of the numbers of moments 
employed-of the two-rebit probability density pi(0). At the final (12,700th) iteration recorded, 

2 

the estimate is 503.17735, while at the 12,600th iteration it was slightly less, 503.17064. 

Estimate of /?i(0) 
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FIG. 5: Estimates-as a function of the numbers (greater than 10,000) of moments employed-of 
the two-qubit probability density Pi(0). At the final (11,445th) iteration recorded, the estimate is 
389.99611337, while at the 11,440th iteration it was slightly less, 389.99611028. 

gence properties in these cases (a = |, 1) than the Legendre-polynomial-based procedure we 
have so far employed. However, our attempts to pursue such a course seemed to yield rather 
unstable estimates. The Legendre-based procedure starts with a uniform baseline density, 
while the Jacobi-based method commences with a beta distribution fitted to the first two 
moments.) 

In Fig. 4 , we show the estimates of the two-rebit density pi (0) as a function of the number 
of moments employed in the inversion procedure. In Figs. [5] and [6j we show the estimates of 
the two-qubit density p\(0) and two-quaterbit density ^2(0) as a function of the number of 
moments exceeding 10,000 employed. In Fig. [7j we show the estimates of the generalized 
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TABLE I: y-intercept Estimates p a (0) 



a # moments 



Estimated ^-intercept" 



Candidate Exact Value 



11 

2 



13 
2 



9 

10 
11 
12 
13 
14 
15 
16 



12,700 

11,445 

10,920 

10,500 

8,100 

8,400 

5,900 

7,500 

5,600 

7,500 

5,900 

7,000 

5,300 

9,240 

7,600 

7,000 

7,000 

8,460 

7,770 

8,500 

6,325 

5,500 

6,100 



503.17 (Fig. 
389.99 (Fig. 
296.34924 



4) 
5) 



218.3823524 (Fig. 
157.34353772 



6) 



390 = 2 • 3 • 5 • 13 



111.536231883 (Fig. 7|) 
78.1142995070 
54.202097676095 
37.335822413322 
25.5666434018285 
17.422692059774 
11.82486375905315966 
7.9980614144240790 
5.393779773089919955165 
2.4351194168151571075324 
1.090697637286080985468617 
0.485464865631244541283680427 
0.21497802182219253282721028797 
0.09479712852190406271015841397642 
0.041653632242518789148830886711329662 
0.01824714034728285015352436495808 
0.0079726370509246020901593058511607 
0.003475526958451101365798473125180273 



7425 
34 

11698020 
74347 

7696 
69 



3 2 -5 2 -ll 

2- 17 

2 2 -3 5 -5-29-83 
7- 13- 19-43 

2 4 -13-37 

3- 23 



14423395 5-7 2 -17-3463 
266104 2 3 -29-31-37 



1959591816 
76646425 



1090539700 
92224293 



2 3 -3 2 -7-13-61-4903 
52.37.41.43.47 



2 2 -5 2 -31-61-73-79 
3-7-41-43-47-53 



1111265181870 
206027169929 

5161595531499 
2119647806944 

241329925246960 _ 
221261985904221 — 

238572701926614 
491431448116025 

8495128297888277220 
39516264155200777043 

662143216303197685 _ 
6984844653287163753 



2-3-5-17-29-75136253 
7 3 -47-53-59-61-67 

3 2 - ll-41-97-13109713 
2 5 -53-59-61-67-71-73 

2 4 -5 ; ll-23-29-37-97-109-1051 

3 3 - 59-61-67-71-73-79-83 

2-3 4 -ll 2 -4M777-167051 
5 2 -67-71-73-79-83-89-97 

2 2 -3 5-7-19-11861033-29917361 
" ll 2 -71-73-79-83-89-97-101-103-107 

5-7 2 -29-61-6719-227380583 
3 2 -79-83-89-97-101-103-107-109-113 



96610235322575857186 
2319371207775705136339 

5306622762208739719059 
290819419438451333002342 

2771311717360093058282 
347602894708307171124675 

6825762937391268808429713 
1963950508510291938204587648 



a Estimates are reported to the number of decimal places that remain unchanged even after the incorporation 
of the last one hundred moments employed. For a < 11, the estimates converge monotonically upward- 
providing lower bounds on true values-an observation we take into account in advancing the candidate 
values. 
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Estimate of p 2 (0) 
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FIG. 6: Estimates-as a function of the numbers (greater than 10,000) of moments employed-of the 
two-quaterbit probability density P2(0). At the final (10,500th) iteration recorded, the estimate is 
218.382352474952, while at the 10,499i/i iteration it was slightly less, 218.382352474705. We note 
that, quite closely, ^ff *j 218.38235294. 

Estimate of pi(0) 
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FIG. 7: Estimates-as a function of the numbers of moments employed-of the generalized 
two-qubit probability density P3(0). At the final (8,400th) iteration recorded, the estimate is 
111.5362318835796, while at the 8,300th iteration it was slightly less, 111.5362318835462. Note 
that, remarkably closely, ^ ~ m.5362318840579. 



two-qubit density p%(0) as a function of the number of moments employed in the inversion 
procedure, while in Fig. |8j we show the corresponding estimates of the first derivative of the 
density ^(O). (In an effort to discern an exact value for at least one a-specific value of p' a (0), 
for a = 4, due to the apparent speedy convergence properties there, we employed 11,480 
moments in the Provost procedure [11] and obtained an estimate of -204971.49610681565 
[continually decreasing in absolute value]. The WolframAlpha website suggested a num- 
ber of possible exact values, such as - V41465479089 ^ « -204971.496106763.) Of course, if 
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Estimate of pi (0) 
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FIG. 8: Estimates-as a function of the numbers of moments employed-of the first derivative p 3 (0) 
of the generalized two-qubit probability density for a = 3. At the final (8,400f/i) iteration recorded, 
the estimate is -289597.1343229, while at the 8,300 tfi iteration it was slightly less, -289597.1343453. 

one had in addition to the exact value of p a (0), derivatives of higher order of p a i\p PT \) at 
|p PT | = 0, one could attempt a Taylor series expansion about that distinguished point. How- 
ever, in practical terms, rates of convergence of estimates-using the Legendre-polynomial 
reconstruction procedure [TT] -become increasingly slow as the order of the derivative in- 
creases (cf. [12J). (Therefore, at the present stage of analysis, we are unable to advance any 
possible, plausible candidates for exact values for these derivatives.) 

In Table [nj we list some estimates we have obtained of the first derivative, p' a (0), for 
various values of a. We do note that these values are all negative (decreasingly so as a 
increases), with the exception of the generic 9-dimensional two-rebit and 15-dimensional 
two-qubit cases, a = |, 1, respectively. Further, since the Legendre-polynomial probability 
distribution reconstruction algorithm was relatively fast in the two-quaterbit (a = 2) case, 
we pursued some additional analyses there, complementing our results of for ^(0) 
(Table |l| and ^ for the cumulative ("separability") probability over [0, ^j. Thus, based 
on 5,000 moments, we obtained estimates of P2(— ^g) and ^(^j) of 83.56342 and 0.020157, 
respectively. Also, our estimates of the cumulative probability over the intervals [— ^s] 
and [5^, were 0.81418669739 and 0.000002538812, respectively. In Fig.gwe attempt a 
companion reconstruction of the (sharply-peaked) two-quaterbit determinantal probability 
distribution P2(|p PT ||)- While the separability probabilities for /ia(/-integral a are, it seem 
abundantly clear, rational- valued, along with the ones for integral a, we presently lack any 
evidence/argument that the values of p a (0) are rational- valued for half-integral a (if for no 
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TABLE II: Estimates of the first derivatives p' a (0) 
a # moments Estimated first derivative 
\ 12,700 6410356. 



10,000 339196. 

8,700 -162646. 

10,000 -289800. 

8,100 -311862.10 



8) 



3 8,400 -289597.134 (Fig 
I 5,800 -249449.43420 

4 11,480 -204971.49610681 
I 5,600 -162971.8709833 

5 7,500 -126442.27058586704 
^ 5,900 -96249.9203477421 

6 7,000 -72156.15107106951 
f 5,300 -53418.9817136894241 

7 5,300 -39134.0668213575364 

8 7,600 -20472.80618300754317070 

9 4,800 -10425.4751799502175604 

10 4,800 -5199.062697559066885103 

11 3,500 -2549.69007308405673430 

12 5,300 -1233.4041330231218163946548 

13 9,500 -589.88846332521988518644909 

14 4,900 -279.4104922820162161974541786 

15 5,500 -131.2569443034131032177602791621 

"Digits to the right of the decimal point are reported to the extent they have remained unchanged after 
the incorporation of the last one hundred moments. If no such digits are given, simply the integral part of 
the last estimate obtained is presented. 
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FIG. 9: The two-quaterbit (a = 2) determinantal probability distribution, based on the first 5,000 
moments Ml) of \p PT \ 




FIG. 10: Joint plot of the logs of the a-specific separability probabilities (lower line) together with 
the logs of the estimated values (upper line), based on 10,700 moments, of the y-intercepts p a (0) 
for the 210 values a = |, 1, §, 2, . . . 105. 



other reason than that we have been predominantly concerned with the integral cases). 



Fig. [10] is a joint plot of the logs of the a-specific separability probabilities together with 
the logs of the estimated values (based on 10,100 moments) of the y-intercepts p a (0) for the 



210 values a 



21 i 2> 



2, . . . 105. The correlation between the pair of 210 log-values is almost 



perfect, that is, 0.99988 (cf. Fig. 11). 
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FIG. 11: Ratio of logs of estimated values of p a (0) to logs of separability probabilities 

III. FISHER INFORMATION EQUALITY CONJECTURE 

The probability distributions p a (\p PT \), discussed above, form a one-parameter family- 
parameterized, obviously, by a. Therefore, we can inquire as to the Fisher information [15] of 
this family. In Fig. [12J we plot our estimates of this quantity using the Legendre-polynomial 
probability distribution reconstruction procedure [TTj based on the first 245 Hilbert- Schmidt 
a-specific moments Q of |p PT |. One can, thus, conceive of a scenario in which one obtains 
a stream of determinants of random (with respect to Hilbert-Schmidt measure) partially 



transposed 4x4 density matrices (|p PT |). Then, employing the Fisher information (Fig. 12) 
as a prior distribution over the Dyson-index-like parameter a, one can derive-making use 
of Bayes' Theorem in conjunction with formulas for p a (\p PT \)-a. posterior distribution over 
a (cf. [IS]). (The extensive, diverse analyses reported in the treatise "Science from Fisher 
Information" [15] are concerned primarily with families of translation-invariant probability 



distributions- which clearly the family p a {\p FT \) is not.) In Fig. 13 we show a plot comparable 
to Fig. [lj but now not employing the Hilbert-Schmidt moments of |p PT |, given by Q, but 
rather the Hilbert-Schmidt moments of \p\ [T71 sec. Ill, Fig. 3] (using Pochhammer symbol 
notation) 

n 256-"(l) n (a + l) n (2a + l) w 
( ' Pl) (3a + i) n (3a + §) n (3a+I); W 
determining a one-parameter family of probability distributions ? a (|p|) over \p\ e [0, ^g]. 
(Dunkl has explicitly constructed the family q a {\p PT \ [191 Ex. 5.2].) In Fig. 14, we show 



the counterpart to the Fisher information plot Fig. 12 Since these two plots appear to be 



quite comparable, it readily suggests the conjecture that the Fisher information is identical 
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FIG. 12: Estimates-based on the first 245 Hilbert-Schmidt moments of |p |— of the Fisher infor- 
mation of the a-parameterized family of probability-distributions p a (\p PT \) 




FIG. 13: The a-parameterized family of probability distributions g Q (|p|) over \p\ G [0, ^gg] estimated 
on the basis of the first 100 Hilbert-Schmidt a-specific moments of \p\ 

for the two families of probability distributions p a {\p PT \) an d q a {\p\)- (The correlation 
coefficient between the 32 values-a = |, |, |, 1, . . . 8-used in the two plots is remarkably 
high, 0.994524.) Since the "balanced" set of moments (JsJ) are in some sense "intermediate" 
between the moments of these two probability distributions (Q, ([6])), one might also further 
speculate that the associated Fisher information is identical to the hypothesized other two 
as well. 

If we turn from the Hilbert-Schmidt measure [H E] to the Bures measure [HI HE] over the 



13 



Fisher information estimate 
0.7 : 
0.6 : 
0.5 : 
0.4 : 
0.3 : 

0.2 ; 

0.1 : 



FIG. 14: Estimates-based on the first 100 Hilbert-Schmidt moments of |p| — of the Fisher information 
of the a-parameterized family of probability-distributions <7 a (|/9|). The correlation with the results 



in Fig. 12 is high (> 0.99). 



generalized two-qubit states, then we have the set of moments (cf. [6]) 

= T(6a + 2 )2- to - 8 "- 1 r (n + |) T (n + a + |) T(2n + a + 1) 
UPI ;Bures yv?r(n + 2« + l)r(n + 3a + l)r (2n + 3« + |) ' U 

Our (Mathematica) estimate of the Hilbert-Schmidt Fisher information for the family g a (|p|) 
at a = | was 1.6595, while, in an independent Maple calculation, C. Dunkl obtained 1.660 
(cf. [19J). Dunkl also estimated the Bures counterpart of this value to be 0.9609. Fur- 
ther, for a = 1, Dunkl obtained 0.32532 for the Bures information, while our Mathematica 
Hilbert-Schmidt result was 0.590666. So, it would certainly seem, as might have been antic- 
ipated, that the generalized (a-parameterizedO two-qubit Hilbert-Schmidt and Bures Fisher 
information are different in nature. 
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